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ABSTRACT 

The oxygen absorption line imprinted in the scattered light from the Earth-like planets has been 
considered the most promising metabolic biomarker of the exo-life. We examine the feasibility of 
the detection of the 1.27 /am oxygen band from habitable exoplanets, in particular, around late- 
type stars observed with a future instrument on a 30 m class ground-based telescope. We analyzed 
the night airglow around 1.27 jam with IRCS/echelle spectrometer on Subaru and found that the 
strong telluric emission from atmospheric oxygen molecules declines by an order of magnitude by 
midnight. By compiling nearby star catalogs combined with the sky background model, we estimate 
the detectability of the oxygen absorption band from an Earth twin, if it exists, around nearby stars. 
We find that the most dominant source of photon noise for the oxygen 1.27 jam band detection comes 
from the night airglow if the contribution of the stellar PSF halo is suppressed enough to detect the 
planet. We conclude that the future detectors for which the detection contrast is limited by photon 
noise can detect the oxygen 1.27 jam absorption band of the Earth twins for ~ 50 candidates of the late 
type star. This paper demonstrates the importance of deploying small inner working angle efficient 
coronagraph and extreme adaptive optics on extremely large telescopes, and clearly shows that doing 
so will enable study of potentially habitable planets. 

Subject headings: astrobiology - Earth - scattering - techniques: spectroscopic 



1. INTRODUCTION 

Oxygenic photosynthesis is currently the predominant 
source of energy for life on Earth: 

6C0 2 + 12H 2 C 6 H 12 6 + 6H 2 + 60 2 . (1) 

To be more specific, oxygen molecules are generated by 
the transfer of electrons from water: 



2H 2 + 4n 7 -> 2 + 4H+ + 4e" 



(2) 



where 7 indicates a photon and n specifies the n-photon 
process of the photosynthesis. Terrestrial oxygenic pho- 
tosynthesis has the n = 2 photon process to utilize the 
main band frequency of the sunlight A = 0.4 — 0.7/im 
(Photosynthetically Active Radiation; PAR). The abi- 
otic production of oxygen molecules is difficult since a 
main source of abiotic oxygen molecules on the Earth, 
the direct photolysis of water molecules, is o nly effective 
for e nergetic photons in UV, A < 180 nm (|Leger et al.l 
l201ll ). Hence oxygen and ozone absorptions have been 
considered as the most promising biosigna tures for the 
search for life on terrestrial exoplanets ([Owed Il980t 
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Ang el et al.l 119861: ILeger et al.l 11993b ISelsis et al.l 120021 : 
iKalte negger et al.l 120101 : Leger et al. 1201 lL and refer- 
ences therein). 

Space-based coronagraph missions proposed so far 
have been aimed at detecting the oxygen 0.76 jam and 
0.69 /am bands from habit able planets main ly around G- 
type stars (e.g. TPF-C iLevine et al.ll2009h . Since it is 
difficult for space missions to utilize large telescope (a < 4 
m aperture is typically assumed) , direct imaging of hab- 
itable planets around late-type stars is more challenging 
due to their small angular separation from the host star. 
On the other hand, 30 m class or larger ground-based 
telescopes have also been under consideration for next 
generation missions including the Thirty Meter Telescope 
(TMT), the European Extremely Large Telescope (E- 
ELT) and Giant Magellan Telescope (GMT). Several de- 
tectors for direct imaging of exo planets have been pro- 
posed for these te lescopes (e.g. Mac intosh et al.l 2006; 
iKasper et al.l l2Q10h . The advantages and disadvantages 
of the direct imaging of the groun d-based telescope have 
been wid ely di scussed so far (e.g. iMountain et al.l 120091 : 
lAnthonvl l200l lHanot et al.ll2010D . While several au- 
thors pointed out that the Extreme Adaptive Optics 
(ExAO) has difficulty to achieve the contrast 10 -8 , which 
is not sufficient for the Earth-like planet detection (e.g. 
Stapelfeldt 2006), several authors claimed that the post- 
processing with ExAO can overcome the 10 -8 limit and 
achieve the contrast below 10 -9 . Moreover, the contrast 
of a habitable Earth-size planet around late-type star is 
much more favorable than that of G-type star (10 -10 ), 
for instance, ~ 10 -7 — 10 -9 for late-type stars since the 
reflected-light luminosity is more or less constant among 
all habitable planets due to the temperature constraint, 
while the luminosity of stars differs by several orders of 
magnitude between spectral types. 
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For ground-based direct imaging, near-infrared (NIR) 
observation is preferable over visible thanks to bet- 
ter ExAO performance. Besides absorptions at visible 
bands, oxygen molecules also have the absorption fea- 
ture in the NIR band, at 1.27 /im with a 0.02 /im width, 
which is the second strongest absorption feature next to 
0.76 /im band. While the detectability of this band by 
ground-based telescope has been studied in the context of 
the transm ission spectroscop y of exoplanets around late- 
type stars ([Palle et al.ll201lh . it has not been considered 
in the context of the reflected light in detail. 

Recently several missions to aim to detect the habit- 
able planets around late-type stars have be en pro posed 
(e.g. iMatsuo fc Tamurall2Q10t iKasper et all 12010ft . The 
key concept to detect habitable planets from the ground 
is the combination of a coronagraph an d subsequent 
post processing (e.g. iHinklev et al.l 1201 lft . Although 
the key technology of ExAO, coronagraphs and post- 
processing have not yet demonstrated the performance 
level required to detect habitable planets, many proposed 
techniques are currently under devel opment aiming to- 
ward reaching the photon noise limit (iPerrin et al.ll2006l: 
Sauyage et all l2006[ iGalicher et all 120101: ICornia et al.l 
20101 : iBaudoz et al.l [20101: IVogt et all l201lL references 
therein). For instance, IBaudoz et al.l (|2010l ) presented 
simulation results of the self coherent camera (SCC) as- 
suming E-ELT and showed that the SCC can achieve the 
photon noise limited when the magnitude is below ~ 7 
mag for the H band. In addition, speckle nulling tech- 
nique, which is not a post-processing method, is also a 
promising method to detect a habitable planet around 
a late-type star. Hence it is crucial for astrobiology on 
exoplanets to know what kind of science at the habitable 
planet is possible by ground-based telescopes assuming 
that an idealistic photon-limited detection is reality in 
the near future. 

In this paper, we focus on the 1.27 /im oxygen absorp- 
tion band imprinted in reflection lights as a biomarker 
of oxygenic photosynthesis and examine its detectability 
with the above 30-40 m class telescopes (Extreme Large 
Telescopes; ELTs) with near- future idealistic instruments 
against the speckle noise. 

The rest of the paper is organized as follows. We first 
describe the outline of the absorption detection and clar- 
ify the main source of noise in §2. In particular we es- 
timate the intensity of the sky background by analyzing 
the real data of the sky on Mauna Kea. With several 
nearby star catalogs, we estimate the feasibility of the 
oxygen 1.27 /im detection assuming an Earth- twin at the 
inner edge of habitable zone (IHZ) in §3. In §4, we dis- 
cuss the availability of the stellar radiation to use the 
oxygenic photosynthesis for our sample. Finally we sum- 
marize our results in 85. 



2. 



DIRECT IMAGING OF EARTH-LIKE PLANETS BY 
GROUND-BASED TELESCOPES 



The detectability of the planet itself and the detection 
of absorption lines are different. The former is in many 
cases dominated by the speckle noise from the leakage 
of the main star, as shown in the left panel of Figure [T] 
(we use the term "leakage" to describe the contribution 
of the stellar PSF halo throughout this paper). Thus, 
the signal-to- noise ratio is / p /ctl, where f v is the flux of 
the planet and ox is the standard deviation of the leak- 



age / p . The sky background is expected to be uniform 
on the detector plane, and its non-uniformity and pho- 
ton noise are negligible compared with the speckle noise. 
On the other hand, once the planet has been detected, 
the photon noise of the skybackground / s k y , the planet 
signal, and the leakage is the main source of the statis- 
tical error of the band detection (right panel of Figure 
[T|). In this section, we concentrate on the detectability 
on ground of the 1.27 /im band assuming that planets 
themselves are detectable. We will discuss the relation 
between the detectability of planets and the oxygen band 
with the nearby star catalogue in §3. 



Photon counts 



Photon counts 




Detector Plane 



AA 



Fig. 1. — Schematic view of the planet detection (left) and the 
absorption band detection (right). The photons from the planet, 
the leakage of the main star, and the sky background (airglow) 
are painted by blue, orange, and green, corresponding to their flux 
Jp5 /l 5 and /sky Green and red arrows indicate the range for the 
absorption bands and for the continuum determination. 



2.1. Signal-to-Noise Ratio of Detection of the 
Absorption Band 

Here we summarize our observational strategy to de- 
tect biomarker absorption bands. The flux within the 
solid angle of the aperture AQ where the planet is lo- 
cated is expressed as 

A (A) = T(A)/ P (A) + / sky (A) + T(A)/ L (A), (3) 

where T(A) indicates the atmospheric transmission, 
/ P (A) is the reflection light from the planet, /l(A) is 
the leakage from the main star (speckle noise and halo), 
/sky (A) = S s kyAQ is the sky background within the aper- 
ture due to the night airglow, and S s ^ y (X) is the surface 
brightness of the sky background. 

Simultaneous observation of the planet and the main 
star is essential to eliminate the effect of atmospheric 
transmission variation. This requirem ent is compatible 
with nulling; coronagraph concepts (e.g. lShao et stl.ll2Q04t 
iMurakami fc Baba 201(J). The flux from the observation 
of the main star is expressed as 

/ B (A) = T(A)/*(A) + / sky (A) « T(A)/*(A), (4) 

where /*(A) is the stellar flux. Since / s ky(A) can be si- 
multaneously obtained on the outer side of the detector 
plane, we can estimate the summation of reflectivity and 
the contrast difference 



J(A): 



/a (A) - / sky (A) 



/b(A) 



/p(A) 
/*(A) 



■Cx(A), 



(5) 
(6) 
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C x (A)e 



/l(A) 
A(A) 



(7) 



Here we assume Cx(A) does not contain absorption-like 
feature. We discuss this possibility in £12.51 

Since the photon noise in /a (A) dominates the statis- 
tical noise of J(A), we define the signal to noise ratio of 
the absorption band as 



(S/N) 



abs 



abs 



\/ A/" s ky + Mp + Nh 



(8) 



where A/" s ky 7 A/" p , and A/l are the photo-electrons of the 
sky background, the planet, and the leakage within the 
absorption band width A A at the band center A c . The 
photo-electrons in the absorption band A/" a bs (Fig. PQ 
right) are related to the equivalent width of the band 

Valine as 



Wi 



line 



AA 



W, 



line 



AA 



(9) 



where Wn n e is the equivalent width of the band. 

Equation (jSJ) implicitly assumes that the continuum 
level, as indicated by red arrows in Figure [TJ can be de- 
termined without uncertainty. While this is not the case 
for real observation, the wavelength width available to 
determine continuum (width of red arrows) is 10 times 
larger than the absorption band width (width of green 
arrow) for the 1.27 /im band (see the mock spectra in 
Bottom Panel of Figure [2]) . Hence we ignore the uncer- 
tainty of continuum level. We note that an OH lines sup- 
presser can reduce this uncertainty further. In Equation 
J5J, we also ignore the statistical noises from the denom- 
inator of Equation (|5j) due to the stellar flux /b(A) and 
the offset observation / s k y (A) since /b(A) is much larger 
photon counts and we can sum up the wider region on 
the detector than AQ for / s ky(A). 

2.2. Reflection from the Planet and Strength of the 
Oxygen 1.27 \im Band 

The reflected light from the planet under Lambert dif- 
fuser assumption is expressed as 

/p(A) = ^(A)(^) 2 A(A), (10) 

0(/3) = [sin/3 + (tt - p) cos/3]/tt, (11) 

where A(X) is the Bond albedo, R p is the planetary ra- 
dius, and a is the star-planet distance. The Larmbert 
phase function is characterized by the phase an- 

gle /3 = /(star — planet — observer). We estimate the 
photo-electron from the planet as 

Ae 

he 



^ p «^rps F r(A c ) 7 ^/ p (A c )AA ) 



K = tt ( y ) T oxp ?- T r Q , 



(12) 
(13) 



where D indicates the effective diameter of the telescope, 
, T exp is the exposure time, rx is the throughput of the 
optical system, tq is the quantum efficiency, and h and c 
are the Planck constant and the speed of light. The tpsf 
is the encircled energy, which quantifies what percentage 
the planet flux is within the analyzed region. We assume 
a circular region of radius X/D, that is, AQ = tt(X/D) 2 ). 



Since Fraunhofer diffraction gives 83% photons in this 
range with Strehl ratio=l, we assume tpsf = 0.7 assum- 
ing Strehl ratio ~ 0.9. 

The oxygen 1.27 /im band of the reflection from the 
Earth has been measured by Earthshine o bservation, 
^1.27/im = 2.42 ± 0.30 nm (jPalle et alJl2009h . The 1.27 
/am band observed in Earthshine consists of mixture of 
the main source O2 lines and the minor sources from 
atmospheric dimers (weakly bounded complex of two 
mole cules by Van del W aals force) of O2 • O2 and O2 • N2 
(e.g. iPalle et all [2009( 1. 

We also estimated the equivalent width by creating 
mock spectra with a radiative transfer code. The plane- 
tary reflection spectra are, in reality, slightly shifted by 
the peculiar velocity of the planet, while the Earthshine 
spectra are not. To take this effect into account, we 
computed the high resolution spectra (dX = 0.01 nm) us- 
ing t he radiative transfer code, libradtran ([Maver et al.l 
2005) with the line-by-line scheme (LBL). We used the 
line-by-line optical properti es for the US-st andard atmo- 
sphere gen erated by genln2 (|Edwardsl ll992) based on HI- 
TRAN96 (|Rothmanlll999h . which is provided by the li- 
bradtran website 9 . We divide the planetary surface to 
192 equal area facets by HEALPix (Gors kiet al.l 12005) 
and compute the contribution from each facet assum- 
ing a planet with /? = 90°. The stellar spectra tem- 
plate for the effective temperature 3500 K and the solar- 
met alkity_J^_j^ed for_the incident flux, which is taken 
from ICoelho et al.l (j2QQ5h . We also assume that d = 5 
pc, 0.4 solar radius for the stellar radius, and the star- 
planet distance 0.15 AU, which corresponds to the IHZ. 
We summed the spectra for the clear sky and cloudy sky 
(optical depth is set to 15) under the ground albedo 0.1 
so as to make the spectra with the cloud cover fraction 
0.5, which yields the albedo of the planet ~ 0.3. 

The planetary spectra are combined with the atmo- 
spheric transmission at the Mauna Kea assuming the pe- 
culiar velocity of the planet v p = km/s and v p = 20 
km/s. The latter corresponds the peculiar velocity of the 
Sun to the mean Galaxy rotation. We consider two ex- 
treme cases for the airmass and the water vapor column 
(WVC) of the telluric absorption T(A), (airmass=1.0, 
WVC=lmm) and (airmass=2.0, WVC=5mm). The 
transmission templates at Mauna Kea site, as shown in 
the top panel of Figure [2j were taken from the Gem- 
ini website 10 , which were generated by the ATRAN 
modelling software (Lord, S.D. 1992, NASA Technical 
Memor. 103957). 

The bottom panel of Figure [2] shows the planetary 
spectra with high resolution (AA = 0.01 thin curve) and 
the low resolution (AA = 2 nm; red) obtained by bin- 
ning the high resolution one. Binned spectra transmitted 
in telluric atmosphere are also plotted for v p = km/s 
(blue) and for v p = 20 km/s (black). The absorption fea- 
ture of the spectra with v p = 20 km/s is slightly stronger 
than that with v p = km/s. Figure [3] enlarges these 
three spectra around 1.27 /im. As shown in Figure [3j 
the peculiar velocity of 20 km/s is larger than the typi- 
cal width of each thin line, which displaces the planetary 
lines from the telluric lines. This is the reason why the 

9 http://www.libradtran.org 

10 http:/ /www. gemini.edu/sciops/telescopes-and- 
sites / observing-condition-constraints /ir-transmission-spectra 
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spectra with the peculiar velocity has slightly stronger 
absorption. However, this displacement of the lines does 
not make significant change to the equivalent width of 
the 1.27 jam band (2.5 nm for v p = Okm/s and 2.7 nm 
for v p = 20 km/s in the range 1.26-1.28 /im) as shown 
in the bottom panel of Figure [2] and both cases have the 
equivalent width ~ 2 — 3 nm, which is consistent with 
the result of th e Earthshine obse rvation (Wi.27/xm = 2.4 
nm, green line iPalle et al.l [2QQ9h . Hence we conserva- 
tively take Wi.27/xm = 2 nm as the fiducial value of the 
equivalent width of the oxygen 1.27 /am band throughout 
this paper. Since the depth of the oxygen bands depends 
on the path of light, mainly driven by clouds, we only 
consider the depth seen in the current Earth and we do 
not consider the other cases. 




1.24 1.26 1.28 
wavelength [//m] 

Fig. 2. — Top panel: Atmospheric transmission (airmass=l) at 
Mauna Kea taken from the Gemini website, which was generated 
by the ATRAN modelling software. Bottom panel: the mock spec- 
tra of the Earth- twin around an M-type star, created by solving 
radiative transfer with the Line-By-Line scheme. The flux with 
the spectral resolution 0.01 nm is shown by a thin line. Thick 
curves are binned in 2 nm interval, corresponding to flux without 
the telluric atmospheric absorption (red), with the absorption and 
zero peculiar velocity (blue), and with the absorption and 20 km/s 
peculiar velocity (black). Figure [3] displays the fine structure of 
these three curves around 1.27 /xm. For referenc e, we also plot th e 
spectra taken from the Earthshine by green line (Palle et al. 2009). 



Although photochemical reactions might alter the 
atmospheric comp osition for different stellar types 
([Segura et al.ll2005h , we ignore these effects since we con- 
sider the aerobic environment and the amount of oxygen 
molecules is unlikely to change in this situation. 

2.3. Night Airglow: Oxygen and OH Emission Lines 
The night airglow is described as 

A4k y ~ K j^(S s ky(\))b An AA, (14) 

where (S s \ /i y(X))b is average of S s ^ y (X) over the band b. 
There are numerous OH lines in the J and H bands, which 
dominates the night airglow. Unless using; a supp resser 
of these lines such as OHS (Iwamuro et all I2001h . the 
average in the J band is (S s ^ y (X))j = fjd\S s k y /A\j ~ 




1.2675 1.2680 1.2685 1.2690 1.2695 

wavelength [/im] 

Fig. 3. — Fine structure of planetary mock spectra around 1.27 
/im in Bottom Panel of Figure [2] The red curve indicates the spec- 
tra without atmospheric absorption. The spectra with atmospheric 
absorption and peculiar velocity of the planet are shown by blue 
(v p = Okm/s) and black curves (y p = 20km/s), respectively. 

10 -15 [J/s/m 2 //im/arcsec 2 ] for the best condition, air- 
mass=l and water vapor =1.0 mm on Mauna Kea 11 , 
where we adopt AAj = 0.17 jam for the band width . 
However, for the sky background around the 1.27 /im, 
one must consider the precise surface brightness around 
1.27 jam with width 0.1-0.2 /im, which includes the atmo- 
spheric O2 emission in addition to the OH lines. Hence 
we analyze t he sky observation u sing archival data by 
IRCS/echelle (|Tokunaga et al .11 19981 ) on Subaru telescope 
in November 24, 2007 (Table CQ). 

We extracted the spectra from 19:12 to 1:34 and from 
3:11 to 5:34 (local time), while we could not use the data 
during 1:34-3:11 since each exposure time in this period 
is too short to detect OH/O2 lines. The spectra were 
obtained with a 0" .56 slit providing a spectral resolution 
of ~5000. This is sufficient for resolving most of OH and 
O2 lines around 1.27 /am. The data were reduced via the 
following procedures. Each frame was dark subtracted 
and flat fielded; the median sky value in photo-electron 
(or the number of photons) along the slit was determined; 
and the spectrum was calibrated in flux (using spectra of 
standard stars) and in wavelength. Bas ed on the OH/O2 
line catalogue of lRousselot et al.l (|2QQ0h . we identify lines 
in the 1.24-1.30 /am range by visual inspection and clas- 
sify them by their origin, the OH line (labeled by red 
color), and the O2 line (blue color) and unresolved mix- 
ture of OH and O2 lines (gray) as shown in Figure SJ 

Figure [5] is the nocturnal evolution of OH (red), 
O2 (blue), and mixture (gray) line flux normal- 
ized at 19:12. While the OH lines show a 
decrease by a factor of 2, which is con sistent 
with pr e vious o bservations (e.g. iRamsav et al.l 
" Content! fl996l: fG lazebrook fc Bland- Hawthorn! 



1992 



2001 



lEllis fc Bland-Ha wthorn 2008), the O2 lines decline by 
an order of magnitude. As a result, the surface bright- 
ness evolution around the 1.27 jam line significantly de- 
creases until the middle of the night as shown in Figure 
[U The significant decrease of the 1.27 jam line strength 
is consistent with the res ult by satellite observation of 
the oxygen 1.27 jam line (|Gao et al.ll201l[ ). Their results 
also indicate that this tendency is more or less univer- 
sal througho ut the year, at l east near the equator (see 
Figure 2 in iGao et ahl 1201 lh . Hence, we suggest that 
the observation aiming to detect the oxygen absorptions 

11 http:/ /www. gemini.edu/sciops/telescopes-and- 

sites / observing-condition-constraints / ir-background-spectra 
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TABLE 1 

IRCS Data for the sky background analysis in Nov 24, 2007. 



Local time (HST) 


ID 




Object 


Airmass 


Exposure / frame 


19:12-1:34 


IRCA00189116 


- 00189182 


SDSSJ01408-0839 


1.1-2.1 


900 s 


3:11-5:34 


IRCA00189237 


- 00189255 


PSSJ1057+4555 


1.2-1.7 


900 s 



5.xl0" i4 
4.x 10~ 14 

r-j 

I 3.x 10~ 14 
$ 2.x 10" 14 



PC 



Cm 



i 



a- 




1250 



1255 
wavelength \pm] 



1270 




OH 

2 
OH+0 2 



1280 1 285 1 290 1295 1.300 

wavelength \pm] 

Fig. 4. — Sky background spectra at Mauna Kea obtained by IRCS/echelle on Subaru telescope. The blue and red curves indicate the 
spectra observed at 19:12 and 5:34. Thin vertical lines specify the line wavelength identified by a visual inspectio n. Blue, red, and gray 
labels beside the vertical lines were identified as the OH, O2, and mixture emission lines based on Rousselot et al. (2000). 

should be taken after 23:00 . We note that the sur- 12 
face brightness in the 1.26-1.275 jam range, which avoids 
the strong OH lines between 1.275-1.280 (see Fig. g]) 
is smaller than that of 1.26-1.28 jam. Taking considera- 
tion airmass uncertainty into account, we decide to use 
the fiducial value of the sky background at the 1.27 /am 
band (5 , sky)i.27±o.oi M m = 5 x 10~ 16 [J/s/m 2 //im/arcsec 2 ] 
for computing A4k y throughout the rest of this paper. 

2.4. Leakage vs Sky Background 

The photo-electron from the sky background and 
the leakage are roughly proportional to S' s ky(A)A^ « 
10- 12 A^[erg/s/cm 2 //im] and Ck(A)/*(A) « C x 
-^Q-5.48-o.4m^ reS p ec tively, for the J-band, where m is the 
apparent magnitude. Hence the sky background domi- 
nates when m > 16.3 + 2.5(p — q), where C — 10 p and 
AQ = 10 g [arcsec 2 ]. For instance, p = —8 and q — —3.5, 
one find m > 5.4, that is, the sky background dominates 
the leakage and one can estimate 

(5/7V) abs ^^^— , (15) 




^Kky + A/p ' 

for most cases. We will discuss the validity of this as- 
sumption using nearby stars catalogs in §3. 

2.5. Possibility of False Positive by Speckle Noise 



Fig. 5. — Nocturnal evolution of the sky emission lines nor- 
malized by value at 19:12 as a function of the local time (Hawaii 
Standard Time). The red and blue curves indicate OH and O2 
lines, respectively. Gray curves are mixture of unresolved lines of 
the OH and O2 lines. The OH line shows a decrease by a factor of 
~ 2 in the OH emission throughout the night, while the O2 lines 
significantly decline by a order of magnitude. 

We revisit the assumption that the starlight speck- 
les do not create an absorption-like feature in equation 
(j6]). The speckle noise appears at the angular separa- 
tion Speckle (A) = X/xo radian from the detector cen- 
ter due to the phase error, where xq is spatial period 
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Fig. 6. — Surface brightness around the 1.27 /im band as a 
function of the local time. The red and blue points indicate the 
surface brightness of the sky background in November 24, 2007 of 
the 1.26-1.28 fim range and the 1.26-1.275 jum range, respectively. 
The latter range is selected so as to avoid strong OH lines present 
in the 1.275-1.280 fim range. The horizontal dashed line indicates 
our fiducial value. 

of a typical sinusoidal ripple of the phase of light (e.g. 
iTraub fc Oppenheimerl l201l[ ). Hence the position of 
speckles depends on the wavelength. We consider how 
the speckle moves in the range of wavelength between 
the oxygen band width A A. Let us describe the distance 
of the planet from the detector center by aX/D in the 
diffraction limit unit, the speckle near the planet is made 
by the ripple of xo « D/a. Then the speckle moves to the 
radial direction by AX/xo = a(AX/ X)(X/ D) in the range 
we considered. Hence, if a < A/ A A the speckle does not 
move beyond the PSF size in the range of AA. For a 
0.1 /im width around the 1.27 jam band ( A/AA ~ 13), 
the speckle does not change for the planet at < 100 mas 
from the center for the 30 m telescope. As shown later, 
the planets detectable by the ground-based telescope are 
generally located closer than 100 mas. In this range of 
planet separation, we can assume that Cx(A) has the 
same shape as the average of contrast Cx'(A). Then 
the wavelength dependence of the contrast should be 
smooth around the 1.27 /im band. The contribution of 
the contrast from the AO smoothly depends on A -2 . Sev- 
eral coronagraphs have smooth wavelength dependence 
around 1.27 /im (e.g. Apodized Pupil Lyot Corona- 
graph; APLC Hinkle v^TaII[20TTl ). or even achromatic- 
ity such as Savart-Plate Lateral-shea ring Interferomet- 
ricNuller for Exoplanets (SPLINE; iMurakami fc Babal 
1201(1 ). With these coronagraphs, we might be able 
to avoid false positives due to speckle noise. However 
speckle patterns also exhibit a degree of chromaticity 
due to the chromaticity of various components of the op- 
tical path that impact the quasi-static speckle pattern. 
Though we ignore the systematics due to chromaticity 
in this paper, we stress again the importance of the sys- 
tematic due to chromaticity. 

3. STATISTICAL ANALYSIS OF NEARBY STARS 

In this section, we examine the feasibility of the 1.27 
/im detection from an Earth twin at the IHZ, if exists, 
around real nearby stars. 



3.1. Nearby Star Catalogs 

We compile several nearby star c atalogs: 
sky catalogue of bright M 



the all- 

Dwarfs Lep ine fc Gaidosl 



20lTh (hereaft er LG11) 

Lep ine fc Sharal l2QQ5h. the 



LSPM catalogue 
Hipparcos and Tycho 



Catalo gs (iPerrvman fcESAl 119971). NStars (lGrave^al.1 
(2003, 2006)), and catal ogs bv lGliese fc Jahreifl (Il991l). 
ICasagrande et all (|2008f ) (C08), and ICasagrande et al.l 
fl2010|) (C10). All data except for LG11 12 were taken 
via VizieR catalogue service ([Ochsenbein et al.l [2000) 
and are compiled so as to combine duplicate entries. 

The stellar flux is computed from the J band magni- 
tude if available (LG11, C08, and C10), otherwise we 
simply use by Planck distribution with estimated stellar 
temperature T* and luminosity L*, 



/*(A) = 



2irhc 2 Rl 
A 5 IP 



exp 



he 



(16) 



where R* = yjL^j (47t<7sb^* ) is the stellar radius, <tsb is 
the Stefan-Boltzmann constant. 

We use T* provided in NStars, C08, and C10, and L* 
for the latter two. For the other catalogs, we estimate 
using the color-temperature relatio n. We use (V-H)- 
T* re lation of the M dwarfs derived by ICasagrande et al.l 
(200 81) for LSPM u p to T* = 4000 K and (V-K)-T* of 
Alon so et al.l (fl996l) (V-K > 0.4) for stars of LG11 and 
LSP M (if not available t he H magnitude), and (B-V)- 
T* by I Alonso et al.l ((1996) assuming th e solar metalicity 
for th e Hipparcos and Tycho Catalogs ([Perrvman fc ESAl 
Il997l ) (B-V > 0.2). In the latter two cases, we restrict 
T* > 3600 K beca use of the approxim ate validated range 
(see Fig la/8a of lAlonso et~atlll996l ). The stellar lumi- 
nosity for T* > 4000 K is computed from the V-band 
magnitude my w ith the bolometric correction ([Flowerl 
[T^lTorresll20Toh . BCy(T^). Since BCy(T*) is not ap- 
plicable for M-dwarfs, we adopt the empirical relation 
of the bolometric mag nitude mhoi and the H-ba nd mag- 
nitude mn derived by ICasagrande et al.l ((2008) to stars 
with T* < 4000 K in LSPM, 



^boi = 1-94 + 1.06 ran- 



(17) 



We assume a telescope on Mauna Kea and eliminate stars 
with declination < —40 degree. Finally we pick up 3294 
stars within 30 pc. The numbers of M, K, G, & F stars 
are 202,36,16,7 within 10 pc and 2299,879,439,262 within 
30 pc, corresponding to the completeness 96, 86, 62, & 
64 % within 10 pc and 41, 78, 63, & 88 % within 30 pc, 
respectively. 

3.2. Planet Assumptions 

For the planet, we assume the current Earth- twin at 
the inner edge of the habitable zone, which has the equiv- 
alent width of the 1.27 fim band Wi.27/im = 2 nm and 
the bond albedo A(X)= 0.3, which is a typical value for 
the Earth. We assume that the planet is at maximum 
elongation f3 = 90 degree. We substitute a by th e inner 
habitable distance given bv iKasting et al.l |l993), 



a = amz 



S e fi 



[AU], 



(18) 



where .S'eff = 1.05 for M and K-type stars 1.41 for G-type 
stars and 1.90 for F-type stars. 

12 http:/ /heasarc. gsfc.nasa.gov/W3Browse/all/md warfasc.html 
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TABLE 2 

Fiducial parameters for observation 



parameter 


value 


planet 


planet radius, R p 


1 Rq 


semi-major axis, amz 


IHZ 


EW of the 1.27 /zm band, Wi. 2 7^m 


2 nm 


the Bond albedo, A(X) 


0.3 


observational condition 


atmospheric transmission, T 


0.95 


airglow at 1.27 /im, (S s ky)l.27±0.0l i^m 


5 x 10~ 16 


( J/s /m 2 //im /arcsec 2 ) 




photon-noise limited detector 


diameter of telescope aperture, D 


30 m 


Inner Working Angle 


X/D = 8.7 mas 


extracted region at the planet, Afl 


2.4 x 10~ 4 


(arcsec 2 ) 




throughput, 7*t 


0.5 


quantum efficiency, tq 


0.75 


exposure time, T ex p 


5 hours for spectroscopy 




1 hour for detection 


raw contrast C ra w 


10 -4 at 10 mas 




10 -6 at 100 mas 


detection contrast 


set by photon noise 




equation (IA1I) 



We use 0.95 as a fiducial value of T, which is typical 
atmospheric transmission at Mauna Kea (0.93-0.96 for 
1.26-1.28 /im) 13 . 

3.3. Instrument performance assumed 

As an ELT, we assume a 30 m telescope at Mauna Kea, 
inspired from TMT. Since the planet detection strongly 
depends on instrument performance, we first examine the 
ratio of the planetary flux vs the sky background. Figure 
[71 shows the ratio 



/p(Ac) 
(5 sky (A)) 6 V2.4x 10 



An 



(19) 



The sky background is brighter than the planetary flux 
for most cases. To detect the planet, the leakage should 
be smaller than the planetary flux since the rms of the 
speckle noise is the same order as the leakage flux. The 
leakage does not affect the photon statistics significantly 
if the leakage is suppressed enough to detect the planet. 
Hence we discuss the feasibility of the absorption detec- 
tion separately from the planet detection ignoring the 
photon counts from the leakage. 

The planet detection is also considered assuming an 
idealistic instrument combined a coronagraph and post- 
processing. We assume the raw contrast C raw = 10 -4 
at 10 mas and C raw = 10 -6 at 100 mas. At least in 
principle, the detection contrast by the postprocessing 
should be set by the photon noise. Since the aim of 
this paper is to show the detectability of oxygen for an 
ideal instrument (as opposed to current instruments), we 
assume a photon-noise limited detection with throughput 
= 0.5. The inner Working Angle (IWA) adopted is the 
diffraction limit. A detailed description of the assumed 
instrument performance is summarized in Appendix A. 

13 http:/ /www. gemini.edu/sciops/telescopes-and- 

sites/observing-condition-constraints/ir-transmission-spectra, 
Lord, S.D. 1992, NASA Technical Memor. 103957 
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FlG. 7. — The ratio of the planet flux and the flux from the sky 
background for the J band ((•S'sky)) within Afl = 2 Ax 10 _4 arcsec 2 , 
which corresponds to the PSF circle for a 30 m telescope. 

For the absorption line detection, we assume a 5 hours 
exposure, which is reasonable from the viewpoint of the 
lowest sky backgrounds in one night as shown in Figure 

m 

3.4. Results 

Figure displays the signal-to- noise ratio for the star 
catalogue computed by equation (j8]) assuming A/l <C 
A/" s ky The linear trend of the S/N - distance relation 
is easily understood: since we are considering the planet 
at the inner edge of the habitable zone, the luminosity 
from the planet is constant and then the S/N is approxi- 
mately proportional to d~ 2 . Since the J-band magnitude 
for the same luminosity and distance depends on the stel- 
lar temperature, the S/N is also a function of the stellar 
temperature, which causes the slight dependence of the 
S/N on the spectral type of the star. While 178 stars, 
which are within 10 pc, have S/N above 5, we stress that 
this estimation assumes that the planet has already been 
detected. 

Here we consider the planet detectability with a pho- 
ton noise limited detection. We theoretically estimate 
the planet detectability for each hypothetical Earth-twin 
with the photon-noise limited instrument for a 1 hour 
J-band observation (see Appendix A and Eq. |A2j ). The 
planets with angular separation smaller than 8.7 mas are 
eliminated due to the IWA. The filled circles in Figure [8] 
show the planets with detection S/N > 5<r in the pho- 
ton noise limit. As written in Appendix A, the detection 
limit depends primarily on the distance and dependence 
on the spectral type is not strong. This distance for 1 
hour is indicated by the vertical dashed line ~ 7 pc. 

In general, the achievable contrast on coronagraphs 
strongly depends on the angular separation, a between 
the planet and the star. The bottom panel of Figure [8] 
shows S/N of the absorption band detection (indicated 
by the size of points) on the angular separation-planet- 
star contrast (c sp = f v / f*) plane. On this plane, we 
display the 5 a detection limit of the planet for proposed 
detectors. Solid curves indicates the detection limit of 
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the photon-noise limited detector for a 1 hour J-band 
observation with TMT. 

Ninety seven stars (79, 17, and 1 for M, K and G stars) 
for the photon-noise limited detector are above the 5 a 
detection limit and almost all these stars (76 stars) have 
S/N > 5 of the oxygen band detection. Since the com- 
pleteness of the M and K stars within 10 pc is high (96% 
and 86%), the number of the stars above the limit in our 
sample is close to the real value. 

For reference, we also plot the 5 a limit for cur- 
rently proposed instruments. The dashed curve corre- 
sponds to the detection limit for the visible-light corona- 
graphic imaging polarimeter IFS (near-infrared integral- 
field spectrograph) for EPICS, though IFS assumes a 42 
m telescope. We also show the 5 a detection limit of 
EPOL (the exoplanet polarimeter for EPICS) at the E- 
ELT for reference although the J-band is not used for 
EPOL. These detec tion limits for EPICS are taken from 
iKasper et all (|2010h . 

In this paper, we do not take single/binary stars into 
consideration since the information about it is quite in- 
complete for the catalog we use. For M and K subdwarfs, 
the estimate d values of the bina ry fraction are around 30 
%: 35 ±5 % dReid fc G izis 1997), 26±3 % (iDelfosse et al.l 
12004( 1. 26 ± 6 % (|Jao et al.ll2009[ ) and references therein. 
Since direct imaging of planets around binary stars is 
difficult, we should multiply our expected number of the 
detectable planets by 0.7. Hence we conclude that plan- 
ets around ~ 70 single stars are above the 5 a detection 
limit and ~ 50 of these have S/N > 5 for the oxygen 
band detection. 

3.5. Planet Detection with the Absorption Feature 

So far, we have discussed the feasibility of the oxy- 
gen absorption assuming the planet is already detected. 
However, it is valuable to consider an absorption feature 
of the planet candidate below 5 a detection since absorp- 
tion feature itself increases the confidence level of the 
planet detection. Statistically, the combination of both 
the 3 a planet detection and the 3 a absorption detec- 
tion provides almost the same probability as a 5 a detec- 
tion. While this situation may be less reliable than the 
5 a planet detect alone, the planet search with spectral 
anomaly such as the absorption feature will be useful for 
large telescopes with good photon statistics. Indeed the 
simultaneous differential imaging (SDI) technique has 
been implemented on several instruments, which utilize 
the 1.6 fim metha ne absorption feat u re for identification 
of gas giants (e.g. iBiller et al.ll2007f ). iBiller et al.l (|2QQ9h 
tested a non-simultaneous spectral differential imaging 
technique (NSDI) which utilize the oxygen 0.76 fim ab- 
sorption, assuming observation of an Earth-like planet 
from space. Since the 1.27 jiva oxygen feature is the 
strongest in J,H, and K bands ex cept for water of the 
Earth reflection (Palle et al. 2009), the utilization of this 
band to detect the Earth-like planet is reasonable. 

4. DISCUSSION 

In this paper, we have considered habitable planets 
with the oxygenic environment created by oxygenic pho- 
tosynthetic organisms. Although the probability of the 
emergence of oxygenic photosynthetic organisms is ex- 
tremely difficult to evaluate, we can compare available 
energy of nearby detectable stars to use the photosyn- 
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Fig. 8. — Detectability of the 1.27 fim band from the Earth twin 
at the inner habitable edge with the 30 m telescope at Mauna Kea 
(declination > -40 degree) and 5 hour exposure. We assume tha t 
the stellar leakage is smaller than the sky background (see fl2.4|) . 
The upper panel displays the signal to noise ratio as a function of 
distance. Filled circles indicate the candidates with detection S/N 
> 5cr for a photon-noise limit (1 hour observation) on TMT and 
with the star-planet angular separation larger than 8.7 mas. The 
vertical line is the approximate detection limit of the photon-noise 
limited detector for 1 hour exposure. The lower panel shows the 
signal-to-noise ratio on the angular separation - the planet-star 
contrast. The 5 a detection limits for the photon-noise limited 
detector, IFS and EPOL for EPICS/E-ELT are shown by solid, 
dashed, and dotted lines (Kasper et al. 20101). 

thesis with that on the Earth. Since the photosynthe- 
sis process on the Earth reduces the photon energy to 
the specific excitation energy at the center of the re- 
action via the light harvesting antenna, the availabil- 
ity of photon energy by the photosynthetic organism is 
quantified by the photon flux, not by the energy flux, 
of the PAR, known as the Photosynthetic Photon Flux 
Density (PPFD). Since the excitation energy at the re- 
action center is ~ 0.7 /im, the PAR is defined in the 
range of 0.4-0.7 /am for the oxygenic photosynthesis on 
the Earth. The PPFD of the habitable planet around 
late- type stars is a half to a tenth of that of the E arth 
re.g. lWolstencroft fc Ravei3l2QQ2t iKiang et al-ll2QQ7h - We 
compute the PPFD for the detectable planets by the 
photon-noise li mited detecto r in ou r sample indicated in 
Figure [8l While IKiang et all ([20071 ) have argued the pos- 
sibility to extend the energy to redder to obtain much 
energy on the planet around the late- type stars, which 
might need to the 3-photon process to obtain enough 
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energy to make oxygen molecules from water, we conser- 
vatively take the same range as terrestrial oxygenic pho- 
tosynthesis (0.4-0.7 jam). Since the central wavelength 
of the V-band 0.55 jam is also the center of the range 
of the PPFD, we extrapolate the V-band magnitude to 
the PAR and compute the PPFD. If the V-band mag- 
nitude is not available, we estimate the PPFD assuming 
the Planck distribution although the visible magnitude is 
systematically smaller than that expec ted from the blac k 
body radiation for the late- type star ([Kiang et al. 2007). 
As shown in Figure the PPFD of most detectable plan- 
ets by the photon-noise limited detector is around 10- 
40 % that of the Earth, which is not significantly low 
for the oxygenic photosynthesis. We adopt the PPFD 
of the Earth from Table A5 of Wol stencroft fc Ravenl 
(2002) with the distance correction (from the IHZ to 1 
AU). The PPFD we computed is the one at the top of 
the atmosphere, not the one at the planetary surface, 
which the organism can use. Cloud coverage and op- 
tical depth are also important to determine the PPFD 
at the planetary surface. These might depend on many 
factors, such as the spin rotation period, water content, 
surface temperature and so on. Observationally, if the 
planet is tidally unlocked, the photome tric variation can 
be used to estimate the cloud coverage (jFord et al. 2001; 
Fuiii et al. 2 0111) or even mapping the cloud distri bution 
(|Kawahara fc Fuiiil f2QlTt iFuiii fc Kawaharal [20121 ). For 
tidally locked planet, the cloud coverage for the illumi- 
nated side might be larger t han that of the Earth (e.g. 
Uoshil 120031 : iKite et al.ll20TTh . Thus, availability of the 
photosynthetic light on planets around the late-type star 
is still an open question for both aspects of cloud for- 
mation and methodology to estimate for tidally locked 
planets. 
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Fig. 9. Photosynthetic Photon Flux Density (PPFD) of the 
detectable planets at the inner habitable zone. The range of the 
Photosynthetically Active Radiation (PAR) is 0.4-0.7 /im. The 
PPFD for filled and symbols are computed from the V-band mag- 
nitude while opened ones are obtained as black body radiation due 
to lack of the V-band magnitude. Circles and rectangles represent 
the planets above 5 a detection by the photon-noise limited detec- 
tor. The horizontal line indic ates the PPFD on the Earth, PPFD0 
(Wolstencroft & Raven 2002). The three dotted lines correspond 
to 40, 20, 10, and 5 % of the Earth PPFD level from top to bottom. 



We did not consider detectability of other important 



biomarkers including H2O, CH4, CO2 in this paper. The 
detection of oxygen with reductant such as methane 
might support the evidence of the existence of exo-life. In 
this case, methane does not have to be biotic molecules. 
It proves that there exists some mechanism which con- 
tinuously produces oxygen, since oxygen is consumed in 
the oxidation of methane. While methan e on the cur- 
rent Earth is more challenging to detect, (Se gura et al.l 
2005) showed that UV photochemical process signifi- 
cantly enhances the absorption of methane for habit- 
able planets around late-type stars ( Figure 8 in their 
paper). As shown in (Seg ura et al.ll2005h . the detectabil- 
ity of methane is complex because the methane is minor 
component in the current Earth (but oxygen is not the 
case). Since the issue that abiotic production of oxygen 
is still under debate, we do not consider the methane 
detection in this paper. The practical problem is that 
we can no longer adopt the simple estimate of S/N since 
these lines are not in the middle of bands and, have wider 
wavelength range. We will consider detectability of these 
biomarkers by generalizing the method in the next paper. 



5. CONCLUSION 

In this paper, the detectability of oxygen absorption, 
as a biomarker from habitable planets, by ground-based 
telescopes has been examined. Analyzing real data of 
the night sky, we showed that the telluric oxygen emis- 
sion as the strong sky background is damped by an order 
of magnitude during a night. With nearby star catalogs, 
which is almost complete for M and K stars within 10 
pc, we found that the oxygen 1.27 jam absorption band is 
detectable with future photon noise limited instruments 
and 30-40 m telescopes for planets around dozens to hun- 
dreds late-type stars assuming an Earth- twin at the inner 
habitable zone around the stars. Photosynthetic pho- 
ton flux density of these planets is not significantly low, 
10%-40% of that on the Earth. We conclude that ex- 
tremely large telescopes in the near future will enable us 
to search for the signature of oxygenic photosynthesis on 
the habitable planets around the late- type stars. This 
paper demonstrates the importance of deploying small 
IWA (~ X/D) efficient coronagraph + ExAO on ELTs, 
and clearly shows that doing so will enable study of po- 
tentially habitable planets. 
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APPENDIX 

THEORETICAL LIMIT OF THE PLANET DETECTION FOR THE PHOTON NOISE LIMITED CASE 

Here we examine the theoretical limit of the planet detect ability for the photon noise limited case. We assume that 
the coronagraph first suppresses the photon of the stellar leakage to N = C raw A/"* at the planet location, where C raw 
is the raw contrast and A/"* is the photo-electron of the main star. 

For the photon- noise limited case, dynamic range is limited by the photo-electron counts of the planet through the 
coronagraph 



dyn = y/rjN, (Al) 

where r] is a detector-dep endent factor. We use rj = 1/2, which is the same value estimated for the Pupil Remapping 
Imager ([Perrin et al.ll2006[ ). Then, the signal-to- noise ratio is expressed as 



(S/N) det = J^ivr = C W ( A2 ) 

The raw contrast is primary determined by the Ex AO. iKorkiakoski fc Verinaudl (|2010f ) performed the simulation of 
ExAO for EPICS and showed that ExAO can achieve the contrast ~ 10 -4 and ~ 10 at the angular separation 10 
mas and 100 mas, respectively (see Fig. 4 of their paper). Hence we assume the raw contrast within IWA as a function 
of the angular separation as 

C raw = max jfc (-^) " 2 , 10- 6 J , (A3) 

where we adopt k = 10 -2 . The condition of the 5 a detection (S/N) det > 5 for the Earth- twin at the IHZ can be 
rewritten as 



2N* fR®Y fL* 



-l 

S eff </>(/3)A(A) > 5. (A4) 



Writing 



w «" w ««(fe)(i4)" ,/(r * ) <A5) 

where Nq^q is the photon count of the Sun observed at the distance of 10 pc with a considered instrument and f (T*) 
is a correction factor due to the stellar temperature dependence between photon counts and luminosity. Equation (|A4[) 
reduces 

d < 2 x lO- 3 [A^ ,io/(T.)] 1/4 pc ~ 7.3 pc, (A6) 

where the last value is derived on the assumption of 5 hours exposure with the instruments in Table [2] and /(T*) = 1.9 
for T* = 3000 K (we assume the black body). 
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